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ABSTRACT 

A numerical  scheme  based  on  polynomial  approximation  is 
presented  for  the  calculation  of  upper  and  lower  bounds  for  the  flexibility 
coefficient  of  end-loaded  orthotropic  cantilever  beams  in  a state  of  plane 
stress,  through  use  of  the  principles  of  minimum  potential  energy  and 
of  minimum  complementary  energy.  The  degr ee-of-freedom  of  the  two- 
dimensional  polynomials  employed  along  with  the  elastic  constants  and 
the  depth-to-span  ratio  of  the  beam  are  external  parameters  of  the 
problem. 
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INTRODUCTION 

A formulation  of  the  problem  of  the  end-loaded  cantilver  beam, 
which  makes  it  possible  to  obtain  upper  and  lower  bounds  for  the  flexibility 
coefficient  through  the  use  of  the  principles  of  minimum  potential  energy 
and  of  minimum  complementary  energy,  has  been  reported  earlier  [l  ]. 

A later  work  [2]  contains  improved  bounds  obtained  through  the  use  of 
arbitrary  functions  of  the  axial  coordinate  in  the  assumed  admissible 
expressions  for  displacements  and  stresses  in  the  calculation  of  the  energy 
functionals.  Although  these  improved  bounds  provide  reliable  corrections 
to  the  elementary  theory  flexibility  coefficient  when  the  parameter  p 
= G/VEE  is  not  too  large,  the  differences  between  the  bounds  becomes 

y 

significant  for  large  values  of  p. 

In  what  follows  we  use  a Rayleigh-Ritz  procedure  in  conjunction  with 
the  Lagrange  multiplier  method  to  improve  the  bounds  obtained  in  [2], 

We  approximate  the  displacement  components  in  the  potential  energy 
calculations  and  the  stress  components  in  the  complementary  energy  calcu- 
lations by  polynomials  in  the  coordinate  variables.  The  coefficients  in 
these  polynomials  are  obtained  by  solving  the  simultaneous  algebraic 
equations  which  result  upon  extremizing  the  energy  functionals  subject  to 
the  constraint  conditions  on  these  polynomials. 
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Numerical  results  obtained  employing  sixth-degree  polynomials  in 
two  variables  show  that  for  isotropic  beams  when  v = 1/3,  and  the 
depth-to-length  ratio  is  as  high  as  0.6  the  maximum  error  in  using  the 
average  value  of  the  bounds  instead  of  the  actual  value  of  the  flexibility 

coefficient  is  less  than  0.  3%.  Using  polynomials  of  the  same  degree  for 

* 

orthotropic  beams  this  error  comes  out  to  be  less  than  1%  when  v = 1/2 
and  p.  = 100  for  all  depth-to-span  ratios. 


FORMULATION 


We  b riefly  restate  the  problem  of  the  end-loaded  orthotropic  beam 
as  formulated  in  [ 1 ].  The  inequality  for  the  flexibility  coefficient  C reads 


* C * 


where  V is  the  prescribed  uniform  end  deflection  and 
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with  2c  and  a representing,  respectively,  the  depth  and  the  length  of 

the  beam,  and  with  E = V EE  . 

m y 

In  equation  (2),  u and  v are  assumed  axial  and  lateral  displace- 
ment  components  which  are  differentiable  functions  satisfying  the  displace- 
ment boundary  conditions 


u (a,  y)  = v (a,  y)  = 0 , v (0,  y)  = - V 


and  in  equation  (3)  CT  and  T are  stress  components  which  satisfy  the 
equilibrium  equations 

CT+T=0,ct+t=0,  (5 

x,  x , y y.y  , x 

and  the  stress  boundary  conditions 

T (x,  ±c)  = 0 (x,  ±c)  = CT  (0,  y)  = 0 (6 

y x 

In  what  follows,  we  obtain  upper  and  lower  bounds  for  C by 

maximizing  I and  by  minimizing  I over  a suitably  restricted  class 
s d 

of  polynomials. 


UPPER  BOUND  CALCULATION 


We  introduce  the  dimens  ionless  quantities 
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X = V G p ' * = ~i ' (9) 
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and  have  therewith,  from  (3),  that  T = V3Ep3I  , where 
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We  satisfy  the  equilibrium  equations  (5)  identically  by  writing 


O = F _ 

x . rjn 
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where  the  A are  arbitrary  constants,  except  for  the  relations 
mn 


E( Zn  - 1 ) A =0  , m — 1,  Z,  . . . , M , 

mn 

n=  1 


y a =o 

mn 


, m — Z , 3,  » • • , M • , 


which  result  upon  satisfying  the  stress  boundary  conditions  r = Cr^  = 0 
at  T)  = ±1,  with  the  condition  a = 0 for  £ = 0 being  satisfied  auto- 
matically. 

Extremization  of  I with  respect  to  the  A , subject  to  the  con- 
s mn 

straint  relations  (13)  and  (14),  is  accomplished  as  follows. 


We  define  a quantity  Ij  by 
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with  B and  C being  Lagrange  multipliers, 
m m 
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Setting  01  =0,  we  obtain 
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where  B,  = 0,  in  view  of  the  form  of  (15).  From  this  we  have  further, 
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and  assume  F in  the  form 
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Equations  (17),  together  with  (13)  and  (14),  represent  a system  of 

equations  of  just  the  right  number  to  determine  the  unknown  constants 

A , B and  C 
mn  m m 

We  next  obtain  an  expression  for  I in  terms  of  the  coefficients 

s max 

A^.  To  do  this  we  make  use  of  the  fact  that  we  may  set  6A^  = in 

equation  (16)  and  write,  using  equations  (13)  and  (14), 

m n aT 

S,  ***.<■“  - 

Considering  the  linear  and  quadratic  terms  present  in  equation  (10), 

we  may  deduce  from  (10)  in  conjunction  with  (18)  that 

.1  _ 

I = r (0,  T))  dT)  . (19) 
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With  T from  (11)  and  (12),  this  becomes 
N 

I = - £ A, 

s max  , In 


We  introduce  the  associated  I on  the  right  of  equation  (1),  and 

s max 

therewith  obtain  the  upper  bound  formula 


£ ”,  A.n 


where  C is  the  elementary-theory  flexibility  coefficient  1 / 2 Ep  . 
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It  remains  to  write  the  system  (17)  explicitly.  To  do  this  we 


introduce  O , T,  and  a from  equations  (11)  and  (12)  into  equation  (10) 
x y 


for  I and  carry  out  the  differentiations  with  respect  to  and  the 

integrations  with  respect  to  £ and  T).  In  this  way  equation  (17)  becomes 
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again  with  = 0. 

Setting  k = 1,  -1=1  in  equation  (22),  we  may  sum  the  series  which 
occurs  in  (2 1 ),  as 


C - 1 
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N 

Aln  = “X* 

n = 1 


(23) 


and  this  enables  us  to  write  equation  (21)  in  the  final  form 


yu 


Xa 


C 1 - C 

o 1 


(24) 


It  remains  to  solve  (22),  (13),  and  (14)  in  order  to  obtain  the  value 
of  the  Lagrange  multipliers  . 

We  note  that  for  M = 1,  N = 2,  the  above  procedure  gives 
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CTT/C  = I + 6X2/5  which  coincides,  as  it  must,  with  a previously 
U o ’ 

obtained  upper  bound  [ 1 ]. 


LOWER  BOUND  CALCULATION 
Introducing  the  nondimensional  displacements 
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Expressions  for  displcements  are  now  assumed  in  the  form 
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where,  since  v (0,  T))  = - 1,  we  have  the  constraint  condition 


X)  amn  + 6ln  0 > n - 1,  2, N . 

m=l 


We  obtain  a and  b by  minimizing  I , with  respect  to  these 
mn  mn  7 d 


coefficients,  subject  to  (28),  that  is  by  extremizing 
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with  the  c^  being  Lagrange  multipliers. 
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Setting  6l^v  = 0 gives 
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Equating  the  coefficient  of  each  6 an<3  6 b i-n  (30)  equal  to 

zero,  we  obtain  the  equations 
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which,  together  with  (28),  allow  the  determination  of  the  unknown  constants 

a , b and  c . 
mn  mn  n 

We  next  obtain  an  expression  for  I,  . first  in  terms  of  the  co- 

d min 

efficients  a and  ultimately  in  terms  of  c alone.  To  do  this,  we  make 
mn  i 

use  of  the  fact  that  we  may  set  in  equation  (30)  6a  . = a . and  6b  = b 
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Since  I is  a quadratic  form,  the  first  two  terms  in  (33)  add  up  to 

21,  . and  we  have  then 

d min 


d min 
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Equation  (34)  gives  with  the  aid  of  the  constraint  equations  (28), 


I , . = c 

d min  l 


which,  in  connection  with  equation  (1),  gives  the  lower  bound 
formula 

CL  , 


It  remains  to  write  equations  (3  1)  and  (32)  in  explicit  form.  We 

accomplish  this  by  introducing  u and  v from  (27)  into  the  equation  (26), 

and  by  carrying  out  the  differentiations  with  respect  to  a and  b 

k-t  k-t 

as  well  as  the  integrations  with  respect  to  £ and  rj  • This  gives 
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Equation  (37),  (38)  and  (28)  must  now  be  solved  in  order  to 

obtain  the  value  of  c in  (36).  Numerical  results  reported  in  the  next 

l 

section  are  obtained  by  setting  M ' = M - 1,  N'  = N,  so  as  to  have  an 

equal  number  of  a-terms  and  b-terms  in  the  experssions  for  the  shear- 

strains  which,  numerical  evidence  indicates,  is  the  most  effective  way 

to  carry  out  this  work.  We  note  that  the  choice  M = 3,  N = ! corresponds 

to  the  "semi-elementary"  polynomial  used  in  the  calculation  of  the  lower 

bound,  C /C  =1  - vs  + Xs  [ 1 ] . 
i— < o 

DISCUSSION 

In  Table  1,  values  of  upper  and  lower  bounds  are  presented  for 
isotropic  beams  when  v = 1/3,  and  c/a  = 0.  3 and  0.  5,  as  functions  of 
the  integers  M and  N representing  the  number  of  degrees-of-freedom 
of  the  associated  admissible  polynomials.  The  Hilbert-type  matrices 
involved  in  the  numerical  work,  require  double  precision  arithmetic  in  the 
solution  procedures  for  large  values  of  M and  N.  The  results  show  that 
in  the  extreme  case  of  a beam  which  is  as  deep  as  it  is  long  ( c/a 
= 0.  5),  the  maximum  percentage  error  in  employing  the  mean  value  of  the 
present  bounds  with  M = N = 6,  instead  of  the  actual  flexibility  coefficient 
is  only  0.  22%. 

Numerical  values  of  the  upper  and  lower  bounds  obtained  for 

orthotropic  beams,  with  the  large  value  u = G/  E =100  when  V 

m 

= 1/2,  are  plotted  against  X = V E/G  c/a  in  Figure  1.  To  obtain  a 
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maximum  error  of  less  than  one  percent  it  is  found  necessary  to  use 


1 


M = N = 6 in  both  upper  and  lower  bound  calculations.  These  results 
are  represented  by  the  curves  U66  and  L66.  For  the  sake  of  comparison, 
we  also  include  the  bound  curves  U2  and  L4  obtained  in  [2].  The 
bounds  U12  and  L31  are  equivalent  to  the  earlier  results  obtained 
in  [ 1 ].  The  bound  curves  U66  and  L66  show  that  there  is  a decrease 
in  flexibility  for  small  X.  As  shown  previously  in  [2],  this  decrease  is 
due  to  the  end  constraint  v = 0 at  x = a.  The  shear  deformation  effects 
dominate  when  X is  large  and  the  flexibility  progressively  increases  with 
increasing  values  of  X. 
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Table  1.  Convergence  of  the  bounds  for 
isotropic  beams  when  v = 1/3. 
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c/a  = 0.5 
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% 
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1.  2805 
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1.  6974 
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2.  19 

4,4 

1.  2260 

1.  2762 

2.  05 

1.  7188 

1.  7632 

1. 29 

5,  3 

1.  2605 

1.  2754 

0.  59 

1.  7447 

1. 7674 

0.  65 

5,4 

1.  2618 

1.  2734 

0.  46 

1.  7457 

1. 761 1 

0.  44 

5.5 

1.  2628 

1.  2732 

0.  32 

1.  7483 

1. 7597 

0.  33 

6.6 

1.  2642 

1. 2715 

0.  29 

1.  750  1 

1 . 7581 

0.  22 

C - C 

1 % Error  = — — X 100. 
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